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Abstract 

We review a self-consistent scheme for modelling trapped weakly-interacting quantum gases at 
temperatures where the condensate coexists with a significant thermal cloud. This method has been 
applied to atomic gases by Zaremba, Nikuni, and Griffin, and is often referred to as ZNG. It describes 
both mean-field-dominated and hydrodynamic regimes, except at very low temperatures or in the regime of 
large fluctuations. Condensate dynamics are described by a dissipative Gross-Pitaevskii equation (or the 
corresponding quantum hydrodynamic equation with a source term), while the non-condensate evolution 
is represented by a quantum Boltzmann equation, which additionally includes collisional processes which 
transfer atoms between these two subsystems. In the mean-field-dominated regime collisions are treated 
perturbatively and the full distribution function is needed to describe the thermal cloud, while in the 
hydrodynamic regime the system is parametrised in terms of a set of local variables. Applications to finite 
temperature induced damping of collective modes and vortices in the mean-field-dominated regime are 
presented. 
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I. INTRODUCTION 



In most experiments involving Bose-Einstein condensates, the system is only partially con- 
densed, meaning that thermal excitations can play an important role in the damping of condensate 
collective modes 03-H) or macroscopic excitations such as solitons [5] and vortices BSEl, making it 
imperative to include the full dynamics of the non-condensed atoms when modelling such systems. 
This includes, beyond the usual mean-field contributions, collisions between the non-condensate 
atoms, and particle-exchanging collisions between the condensate and non-condensate. The ap- 
proach reviewed in this chapter reflects an accurate representation of these combined dynamics, 
and can therefore fully simulate the back-action of the thermal cloud on the condensate, which is 
often neglected, treated to lowest order or in linear response (see Ref. [8J). As this method has 
been implemented by Zaremba, Nikuni and Griffin (following on from early work by Kirkpatrick 
and Dorfman ll9l-fT2"l) we henceforth refer to it as the ZNG method, as described in detail in the 
recent book of these authors lfl~3Tl . 

The strength and frequency of collisions between atoms characterises two distinct dynamical 
regimes IflHl : (i) in the collisionless (or mean-field dominated) regime, in which most experiments 
with ultracold atomic gases are conducted (in stark contrast to helium), the physics tends to be 
dominated by mean-field effects; nonetheless, an accurate description of collisions is essential for 
fully describing the system properties, and accounting for changes in condensate atom number. 
Importantly, a clear separation of timescales (average collision time is longer that the collective 
mode period), enables collisions in this regime to be treated perturbatively (see also Ref. [8]). (ii) 
Some recent experiments have been conducted in the crossover to B3 or deep within the hydro- 
dynamic (or collision dominated) [fT5ll regime. In this regime — which has strong analogies to 
the two-fluid behaviour of superfluid 4 He — , the high gas density leads to very rapid collisions 
between thermal atoms, such that the non-condensate enters a local hydrodynamic equilibrium 
(a precursor to true thermodynamic equilibrium); this enables its description in terms of a few 
local hydrodynamic variables (e.g. local density, velocity, chemical potential, temperature, and 
pressure). 

The power of the ZNG theory lies in its ability (i) to successfully describe both collisionless 
experiments with collective modes lfT6l - PT8Tl and macroscopic excitations |fT9l 1201 at finite tem- 



peratures (see Section IV); and (ii) to reduce to the damped two-fluid equations of 4 He in the 
hydrodynamic limit II2III29H (Section [V]). 
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II. METHODOLOGY 



The second quantised Hamiltonian for a weakly-interacting Bose gas is given by 

r * + IT n 2 v 2 1 ir , * + , , * , ) * 

// = J rfr^(r)| -^ + y -t(r) + -J JrV(rOtf(ry)¥(r')mr), (1) 

where V e xt(r) is the external confining potential and U(r, r') = gc^r - r') describes binary inter- 
actions, with g = Ann 2 aim, where a is the s-wave scattering length. Following Beliaev POll . we 
decompose the second quantised field operator as 

%(r,t) = <t>(r,t) + fi(r,t), (2) 

where we define a condensate wavefunction via the non-equilibrium ensemble average 0(r, t) = 
(^(r, f)), which takes a non-zero value under the assumption of Bose broken symmetry. This in 
turn implies that ($r'(r, t)) = 0, with if/' capturing all fluctuations around the classical mean-field 
of the condensate, thus is often termed the non-condensate operator. 

By taking an average in the Heisenberg equation of motion for the field operator, itid^/dt = 
[*P, H] one obtains the equation of motion for the macroscopic wavefunction (omitting the explicit 
dependence on r and t), 

dd> I h 2 V 2 \ -s + * * 



dt \ 2m 

which is exact in the context of symmetry-breaking. Expanding the final term using ([2]), we ob- 
tain (^F'^FF) = n c <p + m'cf>* + 2n'<p + (ty^yy), where we have defined the following mean-field 
'densities': 

• condensate density: n c (r, t) = \(p(r, t)\ 2 , 

• non-condensate density: n'(r, t) = {if/'\r, f)\j/\v, t)), 

• anomalous pair density: m'(r, t) = (Jj/'(r, t)iff'(r, i)). 

The equation of motion for the condensate therefore reduces to fl3TT - [34]| 

^ + gm>*+g<iA'W'>. (4) 



n 2 v 2 

+ Kxt + g(n c + 2n') 



2m 

Thus the condensate is coupled to higher order correlations which are further coupled to even 
higher order non-condensate correlations, ad infinitum, suggesting the need for a suitable trunca- 
tion. An important feature of this equation, highlighted in 0111331 l34ll . is that condensate growth 
from a zero initial value can only arise through the conventionally-neglected 'anomalous triplet 



term' (tf/^tf/i//), pointing to its importance as a 'source term' in condensate kinetics j|2T|. The 
corresponding equation of motion for non-condensate atoms reads 

h 2 v 2 



2m 



+ Vea + 2g(n c + n) 



(5) 



+ gf (fift - rri) + 2g(j> (frW -ri) + g (<A' W - W>) • 

Note that if all the terms n'(r,t), m'(r,t), and are set to zero, Eq. (H) trivally reduces 

to the zero-temperature Gross-Pitaevskii Equation (GPE). For interacting systems and/or at finite 
temperatures however, these quantities will in general be present and should be appropriately dealt 
with. 



A. Mean-Field Coupling: First Order Effects 

The simplest T > Hartree-Fock mean-field method I135H571 is a limiting case of the above 
equations, in which only normal non-condensate terms involving one creation and one annihilation 
operator are retained, with anomalous contributions neglected in both Eq. Q and the entire second 
line of Eq. ([5]). In this limit, and upon making a semiclassical approximation for the kinetic energy, 
the local energy of the thermal atoms becomes in the Hartree-Fock approximation 

2 2 

§,(r, = + Vext(r) + 2g[n c (r, t) + n'(r, t)] = ^- + C/ eff (r, t), (6) 
2m 2m 

beyond kinetic energy and external potential, this includes a mean-field potential U e n(r, t) = 
^ext( r ) +2g[rc c (r, t) +n'(r, t)] created by the condensate n c (r, t) and non-condensate n'(r, t) densities 
which acts on a thermal atom as it propagates. 

In the Hartree-Fock-Bogoliubov (HFB) extension, all quadratic non-condensate operators are 
maintained in the Hamiltonian, i.e. HFB additionally includes anomalous terms with two like cre- 
ation or annihilation operators; such an approach, which can also be derived variationally ll38Tl . 
is appealing as it relies on a quadratic Hamiltonian, which can be routinely diagonalised by a 
Bogoliubov transformation to a quasiparticle basis Il36l . Despite explicitly accounting for pair 
anomalous averages, and providing a lower total energy for the system, this approach is problem- 
atic as its homogeneous limit leads to a gap in the energy spectrum at low momenta, which violates 
the Goldstone theorem [|39l . This inconsistency can be avoided by neglecting the anomalous av- 
erage altogether (or using other tricks 0O]-|4H), as discussed by Griffin fl32l |45j and implemented 
numerically in Refs. f[46l |47ll . The latter simplified approximation (despite not accounting for 
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many-body corrections to the calculations Il40l l45l 14810 forms a good basis for finite temperature 
perturbative theories, with higher order effects in g essential to account for collisional processes. 



B. Particle-Exchanging Collisions: Second Order Effects 

The important collisional processes describing the transfer of an atom into / out of the conden- 
sate is contained in the triplet correlation {{jj'^ijj'ijj') of Eq. Q I12T1I3T1I341 . Careful consideration 
(see also [|33l|49l) leads to a dissipative Gross-Pitaevskii equation for the condensate wavefunction 
0(r,f) given by EH 

dd>(r,t) ( fi 2 V 2 r , ) 

m at = \~2m~ + Kxt(r) + 8 [ " c(r ' f) + (r ' t)] ~ ^ (r ' 1 0(r ' °' (7) 

where R(r,t) = -(g/rc c )Im(0*(i// t <// i A')X with Im(. . .) denoting the imaginary part, is a non- 
Hermitian source term which allows the normalisation of the condensate wavefunction (f>(r, t) to 
change with time. 

The derivation of the kinetic equation for the non-condensate atoms, based on Eq. ([5]), fol- 
lows from the pioneering work of Kirkpatrick and Dorfman ifTOl [TT| on uniform Bose gases as 
discussed in detail in [fT3l |2T| : here we summarise the main results. Key to this is to impose a 
local semiclassical approximation and to describe the thermal cloud by a distribution function. 
We thus introduce the Wigner distribution function /(p, r, t) for an atom of momentum p, at lo- 
cation r, and time t, defined in terms of the Wigner operator, /, for the non-condensate atoms, 
/(p, r, t) = f dr'e ip r ' /n ^ (r + r'/2, t) ft (r - r'/2, t) , via the expression /(p, r, t) = </(p, r, t)) = 
Tr{p(t,t )f(p,r,t)}; here Tr denotes the trace with respect to the density matrix, p(t,to), whose 
evolution is defined by ifrdp(t,t )/dt = [H eS (t),p(t,t )], where H eS (t) is an effective Hamiltonian 
chosen so as to exactly generate the equation of motion for the non-condensate operator of Eq. ([5]). 
It can be split as H eS (t) = H HF (t) + H'(t), where H HF (t) = f drif/*[-h 2 V 2 /2m + U eS (r,t)]iff' is 
the unperturbed HF Hamiltonian; the perturbation term H'{t) arises from the difference between 
contributions of multiples of non-condensate operators and their mean values; see Refs. [[Rl [1311211 
for more details. The equation of motion for /(p, r, t) thus takes the form 

df(P dt T,t) =^ Tr fa, ^) [/(P, r, to), #eff(0]} , 

= ^Tr fat, t ) [/(p, r, t ), H HF (t)]} (8) 
+ ^Tr{p(t,t ) [f(p,r,to),H'(t)]}. 



The first term in Eq. ([8]) [~ 0(g)] generates the usual free-streaming contributions to the Boltzmann 
equation in the presence of a slowly varying external potential U e s(r, t), with a gradient expansion 
in terms of differentials in position, r and momentum, p; the second term [~ 0(g 2 )] encapsulates 
the particle evolution due to collisions and gives rise to collisional integrals. We thus obtain the 
following kinetic Quantum Boltzmann Equation (QBE) (see e.g. fl2TJ l50l - [52ll ) 

% + ~ ■ V r / - (V r t/ eff ) • (V p /) = CnU, 0] + C22UV, (9) 
ot m 

the right hand side of Eq. Q contains two collisional integrals representing different types of 
collisions between condensate and non-condensate atoms 11914121k here thermal excitations are 
treated semiclassically via Eq. ([6]). 

The collisional integrals involve binary atomic collisions which lead to the scattering of parti- 
cles from initial to final states, where fi = /(p,-, r, t) is the statistical factor for the destruction of a 
particle in state i and (fi + 1), for the stimulated creation of a particle in state i. They are defined 
as follows: 

CnU , = ^ J^J^J + * - P, - p.) 

x d(s c + e 2 -h~ e A )(2nhf[8{9 - p 2 ) - 6(p - p 3 ) - <5(p - p 4 )] 

(10) 

X[(/ 2 + l)/ 3 /4-/2(/3 + l)(/4+l)L 

refers to a collision involving one condensate atom with one non-condensate atom and results in 
atom transfer between the two subsystems, and 

h J (inhy J {inhy J (2nhy ^ 

x 5(s + s 2 -s 3 - s A )[(f + l)(/ 2 + D/3/4 - // 2 (/ 3 + 1)CA + 1)] 

refers to a collision involving the redistribution of two non-condensed atoms. 

In the absence of the collisional terms C n and C 22 in Eq. ([9]), the equation is known as the 
Vlasov equation, arising in diverse fields of physics such as plasma physics, condensed matter 
physics and astrophysics. C22 is the usual collisional term appearing in the ordinary bosonic 
Boltzmann equation (in the absence of condensation). Setting Cn\f\ = 0, as appropriate in 
thermodynamic equilibrium, implies that particles are distributed according to the Bose-Einstein 
distribution — » rc BE (£,) = \e^ e> ~^ - l]" 1 . In the presence of a condensate, C\ 2 = when the con- 
densed and non-condensate subsystems are in local 'diffusive' equilibrium, implying that there is 
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no net transfer of particles between the two subsystems on average. In cases out of equilibrium, 
this term acts to transfer atoms between the condensed and non-condensed parts of the system, 
with overall particle number conservation ensured by the source term in Eq. |7]) being given by, 



h C dp 



^0 = ^777 ^ 7^C 12 [/(p,r,O,0(r,O]. (12) 



The delta functions of Eqs. ( [TO] ) and ( |TTj ), enforce conservation of energy and momentum; 
condensate atoms appearing in C\i have energy s c = + /i c , and momentum mv c , where /j c 

is the condensate chemical potential. Although the exact expression for /i c depends on pair and 
triplet anomalous averages [2T], within ZNG their static values are not included in /u c or in the 
excitation energies (which are both treated to first order in g), thus yielding the reduced expression 

fi 2 V 2 VnJr t) 

/i c (r, t) = — + V ext (r) + gn c (r, t) + 2gn'(r, t). (13) 

2m yn c (r, t) 

It is crucial to note that this does not imply that the simultaneous annihilation or creation of more 
than one particle are ignored, since such terms are implicitly maintained (in an approximate, yet 
consistent manner) to second order in g in the dynamical equations, thus giving rise to the colli- 
sional integral Cn- While this procedure is approximate, the method used here is internally self- 
consistent, and eventually leads to a theory which conserves both total energy and total particle 
number. 

To make contact with physical variables, we note that the non-condensate density can be re- 
constructed from the distribution function via 

n '^= f 7^^P' r ' ?) - (14) 
J (2nfry 

Equations ([7]) and ([9]) are the closed set of ZNG equations for a condensate coexisting in a trap 
with a cloud of thermal excitations. Their solution governs the evolution of the physical variables 
of interest. 



C. Numerical Implementation 

Equation ([7]) is solved in much the same way as the simpler, T = GPE, with the non- 
condensate mean field and the source term calculated from the QBE at each time step; we briefly 
outline the numerical procedure for solving the QBE Q in the mean-field-dominated regime, 
where the full distribution function /(p, r, t) is required to accurately describe the thermal cloud 
— see ll53l for details. 
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We simulate the evolution of the thermal cloud phase space distribution by means of a set of 
test particles. This allows a Lagrangian approach to be adopted, whereby the motion of a phase 
point is followed in time according to Newton's equations of motion dTi(t)/dt = Pi(t)/m and 
Pi(t)/dt = -Vt/ e ff(r, t). These are initially distributed in phase space according to /(p, r, t) - 
y{2nfif Yr=\ _ r ;'(0)£(p _ Pi(t)), where y = N'/N tp is a weighting factor which ensures that the 
physical number N' of thermal atoms is being represented. As test particles are only used in order 
to accurately sample the distribution, there is no need for their number to match the actual number 
of thermal atoms, which can be smaller or larger, with the only requirement that N tp needs to be 
sufficiently large to minimise the effects of using a discrete particle description. 

The density at a particular grid point can be determined by using a cloud-in-cell method ll54l . 
This method takes into account the actual position of a particle within a cell, and also allows for 
migrations between grid cells. As residual fluctuations remain, additional smoothing is necessary, 
achieved through convolving the densities with a Gaussian of width a few times the thermal cloud 
grid spacing (the condensate density is also convolved for consistency except in simulations in- 
volving solitons and vortices, where the main physical effects would become smeared out). At 
the end of each time step, the probability that a given test particle suffers a collision is determined 
and its position and velocity are updated accordingly (respecting momentum and energy conser- 
vation). The collisions are treated using a Monte Carlo sampling technique: whether a collision 
has occurred, and what its type is, is determined by comparison of the calculated collisional prob- 
ability with a random number. Finally, the condensate evolution is appropriately amended by the 
occurrence of any collisions which transfer particles into/out of the condensate. As collisions are 
dealt with separately from the dynamical evolution, we are able to choose which collisional effects 
to include in our simulations, and thus assess the importance of each contribution in a physical 
setting, thereby enabling a comparison to other theories (in appropriate limits). 

D. ZNG Condensate Hydrodynamics 

Introducing phase and amplitude variables to the condensate wavefunction 0(r, t) = 
Vn c (r,?y e(r '° allows Eq. <g to be recast in hydrodynamic form: 

dn c d\ c I 1 ?\ 

— = V • (n c v c ) = -T l2 [f,(/>] and m— = -Vhu c + -mv, I . (15) 
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The condensate velocity is given by v c (r, t) = (h/m)V9(r,t), fi c is defined by Eq. (13) and 
Fn[f,<p] = r 12 (r, t) = 2n c R(r,t)/h is the 'source' term due to particle-exchanging collisions be- 
tween condensate and thermal cloud. 



HI. VALIDITY ISSUES 

A. Validity Domain 

A key feature of a successful theory for describing bosonic quantum fluids is its ability to ex- 
plain the phenomenon of superfluidity. This arises naturally within ZNG, since the condensate 
velocity v c is precisely the velocity describing the flow of the superfluid density in the two-fluid 
model. Although the ZNG equations are based on the Beliaev scenario, this is done here within 
certain approximations for the excitations, namely that they are single-particle-like (Hartree-Fock) 
and they ignore anomalous averages. Thus, the present implementation is inappropriate for de- 
scribing collective phonon-like excitations which become important at very low temperatures. 
This could nonetheless be remedied, in favour of Bogoliubov quasiparticles, since the nature of 
the excitation spectrum is intrinsic to the choice of the Hamiltonian on which perturbation the- 
ory is applied; hence the formalism presented here could also be systematically extended within 
the Beliaev field-theoretic approach to account for such features (see Appendix A of Ref. ll2T1l . 
Chapters 3 and 17 of Ref. O and Ref. (H). 

The ZNG theory has been remarkably successful in describing a range of dynamical problems, 
across an extremely broad temperature range, in both the collisionless and hydrodynamic regimes; 
it even describes the rather involved problem of condensate growth ll55l . provided an initial 'seed' 
condensate is assumed. However, its present implementation will fail in the regime of critical 
fluctuations, i.e. very close to the transition temperature in 3d systems, or over broader temperature 
ranges in ld/2d systems which exhibit enhanced phase fluctuations. 

B. Relevance to Other Theories 

The ZNG theory is formally related to the approach of Walser et al. [|56ll57l . Kohler et al. Il58l 
and Proukakis et al. 0111341 1491 . which amounts to a similar perturbative treatment starting how- 
ever from an appropriately generalised basis which includes the pair anomalous average. However, 
an important advantage of ZNG is its formulation in terms of phase space distribution functions, 
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which allows for a relatively straightforward numerical implementation, that can be extended to 
the hydrodynamic regime. (The closely-related explicitly number-conserving approach ll59l46TTl 
has not yet been formulated in terms of a self-consistently evolving thermal cloud.) 

This book also describes a number of stochastic classical-field approaches (see e.g. Il62l - [65l0 . 
In some sense, the ZNG theory can be thought of as the analogue of the full theory of Stoof ll62l . 
when derived within a symmetry-breaking perspective. The key difference of the ZNG formu- 
lation compared to those of classical field methods is that the ZNG approach directly describes 
the condensate order parameter, whereas the field appearing in the GPE-like equation of classi- 
cal field theories describes a range of 'classical' or 'coherent' modes which include, but are not 
restricted to, the condensate mode. In such theories, the latter is to be extracted (in the Penrose- 
Onsager sense) by numerical diagonalisation a posteriori. In contrast to this, the ZNG method 
directly generates two distinct components, the condensate and the non-condensate (see Fig. [2]), 
and thus it provides an intuitive picture of the fundamental physical processes taking place in a 
partially-condensed system; this also guarantees a direct link to the superfluid properties of the 
system — which are harder to extract from unified stochastic approaches — thus leading directly 
to established theories of superfluid helium (see Section[V]). On the implementation front, the cur- 
rent importance of ZNG is that it includes a fully dynamical thermal cloud self-consistently; this 
is crucial, e.g. in order to correctly predict the Kohn mode induced when suddenly displacing the 
trap, whereas stochastic methods describing a band of modes coupled to a static thermal cloud 
would actually lead to artificial damping. 

IV. APPLICATIONS 

A. Damping of Condensate Scissors Mode 

By analogy with earlier work in nuclear physics (see [66] for a review), the superfluid nature of 
a BEC can be probed by observing the so-called 'scissors mode'; this can be excited [67J by adia- 
batically rotating an axially symmetric trapping potential of the form V e xt(r) = m(a) 2 ± p 2 + to 2 z 2 )/2, 
through a small angle 9 about the v-axis, and then suddenly in the opposite direction through 
an angle -26 . The observation of a single frequency signals the irrotational flow of the conden- 
sate. This should be contrasted to the superposition of two frequencies for the thermal component, 
whose flow has both rotational and irrotational character. Simulations with the ZNG theory by 
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Jackson and Zaremba IfTTll found excellent agreement for T < 0.8T° with a subsequent exper- 
iment measuring the temperature dependence of the frequency and damping of such oscil- 
lations. The inclusion of collisional processes associated with both Cn and C 2 i is crucial for 
explaining the experimental observations, and simulations performed without collisions result in 
up to 50% lower damping rates. The Landau damping process is thus a product of both collisional 
and mean-field effects. 



3 



3.0 

2.8 

2.6 H 

2.4 
0.3 

0.2 

0. 1 

0.0 



o o offii gug I 



* $ o o ^ 



1 T 




i i i i r 
0.0 0.2 0.4 0.6 0.8 1.0 1.2 

T/T c ° 



FIG. 1: (a) Frequency and (b) damping rate of the scissors mode for a variable total number of atoms, 
intended to simulate the experiment of Ref. [3]. The condensate mode is indicated by open (theory) and 
filled (experiment) circles. The open squares in (b) show the calculated average damping rate of the two 
thermal cloud modes, while the filled squares are the corresponding experimental values. Reprinted with 
permission from B. Jackson and E. Zaremba, Finite-temperature simulations of the scissors mode in Bose— 
Einstein condensed gases, Phys. Rev. Lett. 87, 100404 (2001) H3. Copyright (2001) by the American 
Physical Society. 
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B. Decaying Vortex Dynamics 

The thermal cloud leads to dissipation of an off-centred vortex at finite temperatures, causing 
it to minimise its energy in a harmonically-trapped condensate by moving out of the condensate 
radially (6113. Figure [2] (right) shows the nonlinear increase of the decay rate with increasing 
temperature, clearly highlighting the role of all collisional processes, particularly at higher tem- 
peratures Il20ll . The largest increase in damping rates arises from the particle-exchanging Cn 
collisions. Our results are compared to the analytic predictions of Fedichev and Shlyapnikov ll68l 
(FS), which account for the scattering of thermal excitations from the vortex core, and those of 




0.2 0.4 0.6 

T/T 

c 

FIG. 2: Left: Density cross sections of the condensate (top) and thermal cloud (bottom) (T = 0.57c) 
showing clearly the two distinct components simulated within ZNG. Colours range from black (low density) 
to white (high density) with different scales for the two images. Right: (Inset) Log-linear plot of vortex 
radial position as a function of time for T = 0.77c (solid line); the dashed line is an exponential fit, r v {t) = 
r^e 7 ', to the data over < oj±t < 50. (Main Plot) Values of y based on different levels of approximation: 
static thermal cloud with condensate dissipation, iR, included (solid circles), thermal cloud allowed to 
evolve within the QBE (Eq. ([9])) without collisions (open circles), with only thermal-thermal collisions 
(pluses), with only particle-transferring collisions (stars) or with all collisional processes (squares). For 
comparison, analytic predictions of FS ll68l (solid line) and DLS ll69l (dashed line) are shown. (Calculations 
performed at fixed condensate number, N c 10,000, and an initial radial vortex offset ro - 0.26/?tf from 
the trap centre). 
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Duine, Leurs, and Stoof [|69l (DLS), which include the effects of Cu collisions within the static 
thermal cloud approximation; the latter are in full agreement with our results, in the corresponding 
limit of condensate dissipation from a static thermal cloud. The enhanced decay predicted by the 
full theory highlights the importance of including all dynamical processes in modelling experi- 
ments. (See also Ref. IPT91 for the corresponding dark soliton analysis accurately reproducing the 
experiment @.) 

V. RELEVANCE TO OTHER SYSTEMS 

By construction, based on its symmetry-breaking formulation, the ZNG approach is consis- 
tent with prevailing theories for the description of superfluid helium (see also Ref. lfT3lO . In the 
hydrodynamic (or collision-dominated) regime, when the thermal cloud enters a state of local 
hydrodynamic equilibrium, Landau's famous 2-fluid equations for the superfluid and normal den- 
sities arise as a special case of the ZNG equations |[T3ll2Tll22l| . when the condensate and thermal 
cloud are in diffusive local equilibrium with the same chemical potential. These equations can be 
extended to include hydrodynamic damping, which arises when the distribution functions deviates 
from its local equilibrium form, giving rise to the Landau-Khalatnikov two-fluid hydrodynamic 
equations [|23l , of liquid helium, with the correct transport coefficients, thus facilitating a precise 
determination of the crossover between the mean-field collisionless and two-fluid hydrodynamic 
regime. It would be of interest to extend the ZNG formalism to other situations such as dipolar 
Bose gases, spinor condensates, bosonic mixtures and apply this scheme to trapped superfluid 
Fermi gases. 
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